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Abstract 

We prove that the weak version of the SPDE problem 



dV t (x) = 

MO) = o 



-M{x) + -{* 2 M + o*)v;'{x) 



dt - a M Vt (x)dW t , x>0 



with initial density Vq a bounded function and W a standard Brownian motion, has a unique 
solution in the class of finite measure-valued processes. The solution admits a smooth den- 
sity process which has a probabilistic representation and some degeneracy near the absorbing 
boundary. We describe the precise order of integrability of the density and its derivatives near 
the origin via the behaviour of a two-dimensional Brownian motion in a wedge of specific angle. 
Our results are sharp in the sense that we demonstrate better regularity is unattainable. 



1 Introduction 

Let W be a standard Brownian motion on a complete filtered probability space (SI, J 7 , P) and Vq 
be a deterministic bounded probability density function on the half-line (0, oo). In this paper we 
consider the problem of finding a finite-measure valued process v = {^t)te[0T] satisfying the weak 
stochastic partial differential equation 



(4>, 



, v t ) = (4>, Vq) + J fi {<f/, v a ) + X - (a 2 M + aj) (<f>", v s ) ds + a M J (<f>', v s ) dW t 



(1) 



for all t G [0, T] and (j) G C test . We write (0, vt) = f °° (j) (x) ut (dx) and denote the space of bounded, 
twice- different iable functions which vanish at the origin and have bounded derivatives by (7 test . The 
coefficients /i G R, om > and oj > and the time horizon T > are all deterministic constants. 
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Stochastic partial differential equation (SPDE) problems of this type with zero boundary condi- 
tions have been studied extensively |7j , |8] , [10] . The simplicity of our setting allows greater insight 
into the behaviour of the solution near the absorbing boundary. We take advantage of a straight- 
forward smoothing technique to present a more elementary study of the problem than is currently 
offered in the standard literature. We are able to prove uniqueness of solution to ([I]) in the broad 
class of finite-measure valued processes. When referring to "uniqueness" throughout this paper, 
we shall always mean that if v 1 and u 2 are two finite-measure valued processes solving ([!]), then 
E [ {((f>, u t ) - (</>, v t )f dt] = for every <f> £ C7 test . 

The results closest to those presented here are found in [8]. In this paper Krylov demonstrates 
the existence of a unique solution to ([I]) in the class of twice-differentiable function- valued processes, 
u = (tif {x))t£\QT\i w hidi satisfy 



E 



I II 2 , 



I 'II 2 1 
NI2 + 



\xu+ 



dt 



< 00. 



Furthermore if the initial density is iV-times differentiable with x h-> x k V^ E I? (0, 00), for k 
0,1, ... ,N, then it is shown that this solution is (N + l)-times differentiable with 



E 


if 


x k ~ l uf ] 


2 

dt 




Jo 




2 



< 00 



for Jfe = 0,l,...,JV+l. 

A key ingredient in [8], and the majority of the literature on this problem, is the introduction of 
weighted Sobolev spaces. These spaces enable one to control the rate at which the solution and its 
derivatives are allowed to blow-up near the origin. In our one-dimensional setting these weights are 
easy to visualise and it suffices to use the functions 

w c (x) := x c e~ x 

for x G (0, 00) and c£R, The main result of this paper is to investigate how small one can take c 
in the above weighting function whilst preserving square-integrability. The critical quantity is the 
parameter 



7T 



a r- 



+ arcsin 



(2) 



which is connected with the underlying dynamics in the probabilistic construction of the solution to 
([TJ presented in section 2. One should observe that the quantity ix ja — 1 appearing in the following 
theorem is in the range (0, 1). 



Theorem 1.1 (The Regularity and Uniqueness Theorem) Let N > 0. Suppose that Vq is 
bounded and N-times weakly differentiable in (0, 00) with w^p^Vo, ^1-/3/2^0'; w 2—p/2^o ' > • ■ •> w N-j3/2^Q N ^ 



2 



all in 1? (0, oo) for every f3 G (—00, n/a — 1). Then there exists a unique finite measure-valued pro- 
cess, v, satisfying |Ip. 

Furthermore for almost all (co,t) G O x [0, T], vt has a density Vt that is (N + l)-times weakly 
differ entiable with 



E 



T 



W k -l-p/2 V t 



(k) 



dt 



< 00 



for k = 0, 1, 2, . . . , N + 1 and every j3 G (— 00, 7r/a — 1). 



Theorem 1.1 shows that we gain an extra order of integr ability near zero and differentiability in 
the space dimension as the system evolves forward in time. If one takes &m = then ([I]) becomes the 
deterministic heat equation with zero boundary condition. With non-zero gm the solution behaves 
as a heat decay with random fluctuations in location which are driven by the Brownian motion W 
- see Figure 1. This heat flow behaviour is responsible for the smoothness properties of the density 
away from the boundary. It is shown in [9], Lemma 3.1 that the density of the solution is infinitely 
differentiable in the space variable. We can deduce integrability properties of the solution using a 



comparison with the heat kernel (Proposition 2.5). 



As the initial density is bounded and integrable, requiring w^—p^Vo G L 2 (0, 00) is not a 
restrictive condition and allows the derivatives of Vo to be singular at the origin. Notice that since 
Vb is bounded and (3 < 1, the condition w_p/2Vo £ £ 2 (0, 00) is always satisfied. Requiring the 
condition on Vo holds for arbitrarily negative values of f3 does not substantially reduce the class of 
admissible initial densities, since the weighting functions decay exponentially at infinity. We therefore 
allow initial densities whose derivatives grow sub-exponentially in magnitude towards infinity. In 
Theorem |1.1| the case (3 = corresponds to the aforementioned result in jS], and hence our results 
expand upon the known regularity of the solution at the boundary. 

Naturally one would like to know whether we can take /3 > n/a — 1 in Theorem |1.1| It is shown 
in [6], Lemma 5.1 that the critical estimate in Theorem 2.6 required in the proof of Theorem |1.1| is 
sharp, and in turn this allows us to show that the regularity is also sharp: 

Theorem 1.2 (Converse result) Assume that Vq G L 2 (0, 00) is bounded and N -times weakly 
differentiable in (0, 00) with w_pi /2^o, Wi^^^Vq, ^2-/372^0"' ■ • •> w N-/3'/ 2^0^ a ^ * n L 2 (0, 00) for 
every f3' G (— oo,n/a — 1). (One can allow Vq to be smoother by extending this range of (3' to R and 
the conclusion is unchanged.) Then for every /3 G {ir/a — 1, 00) 



E 



T 



w k-l-f3/2V t 



(k) 



dt 



/orJb = 0,l,...,JV. 



The blow-up of the second derivative is an interesting feature of SPDE problems that one does 
not observe in the deterministic setting - that is when o~m = 0. In [8j, Example 1.2 a system similar 
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to ([!]) is presented for which the solution does not have second derivative that is in L 2 (0, oo). An 
explicit description of the blow-up is given in [9] . Suppose \x = and shift the density of the solution 
by setting 

f(t,x) :=V t (x + a M W t ). 
Then / solves the heat equation (with drift) in the region 

{(t,x):te [0,T],x> -<T M W t } 

with zero boundary condition. It is shown (Theorem 5.1, [9]) that there exists A 6 (0, 1) such that, 
with probability one, there is a dense subset of [0, T] on which 

lim t — = oo. 

On this set it is clear that the density can have no space derivative at the absorbing boundary. 

Our study of this SPDE problem is motivated by |3], in which a structural model is proposed for 
the loss due to default in a large portfolio of credit risky assets. The authors construct a finite system 
of Brownian motions, correlated through a common diffusion term, which are killed upon hitting 
a boundary at zero. As the number of particles in the system approaches infinity, the empirical 
measure of the particles - that is the proportion of particles falling in a given set - converges to the 
probabilistic solution v of ([I]) defined in Section 2. That ([I]) has a unique measure- valued solution 
allows us to characterise the model in [3] through the SPDE and Dirichlet boundary condition. 

In the proceeding section we present the construction of the solution to ([I]) and demonstrate 
its essential properties. A proof of the significant (3 + /3)-condition (Theorem |2.6[ ) is presented in 
Section 3. The remaining sections are devoted to proving the above two theorems. 

2 The Probabilistic Solution 
2.1 Construction 

Let us introduce another complete probability space (fl 1 , J- 1 , F 1 ) and independent standard Q 1 - 
Brownian motions W 1 and W 2 . We will denote expectation under F 1 by W . Let the processes X 1 
and X 2 be defined as 

X\ := X t + f it + a M W t + aiWl 

where Xq,Xq are independent jV-random variables with common law Vq. We stop the particles 
upon hitting zero by setting 




x\ te[0, n ) 
d te[n,T] 
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Figure 1: Presented on the left is the evolution of the system in the deterministic case o~m = 0. The 
value of the density at a space-time point is represented by a colour increasing from blue (Vt (x) = 0) 
through to red (Vt (x) = 1). The initial condition is a step function. Compare this with a realisation 
of the evolution for some o~m > 0, which is presentated on the right. We observe heat-type dispersion 
but with a non-smooth fluctuation. 

where d is some cemetery state (one can take d = for our purposes) and 

U :=inf{se [0,T] : X\ < o} . 

Clearly X 1 and X 2 are x f^-random variables, and conditional on knowing X 1 (oj, •) 

and X 2 (oj, •) are independent f^-random variables. It is therefore immediate that X 1 and X 2 are 
equal in f^-distribution, so we make the following definition: 

Definition 2.1 (Probabilistic Solution) Let the probabilistic solution, v, be the measure-valued 
process defined by 

u t {S) :=P 7 (X}GS) 

for all t £ [0, T] and S C (0, oo) Lebesgue measurable. Notice that v is a 0,-random variable. 

Throughout this paper "5 C (0, oo) measurable" shall always refer to a Lebesgue measurable subset 
S. The following useful: 

Proposition 2.2 For all t £ [0,T] and S C (0, oo) measurable we have 



E 



E 



E 1 



L x^s,xfes 
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Proof: Since X 1 and X 2 are equal in f^-distribution we have ft{S) = E 1 
Hence 



L xies 



for i = 1,2. 



E 



MS)' 



E 



E 1 



1 x}es 



E 



E J 



1 xleS 



E 1 



L x?es 



E 



E J 



^xleS,xfeS 



as required, where the final line is due to the f^-independence of X 1 and X 2 . 



□ 



It is a matter of applying Ito's formula to see that v is indeed a solution of Q: 
Proposition 2.3 There exists a full subset of Q on which 



(0,^} = (0,Vb} + 



f fx {(j)', v 8 ) + \ (a 2 M + a 2 ) {4>",v a ) ds + OM [ {4> ,v s ) 
o L 1 J Jo 



dW,. 



test 



for all t G [0, T] and G C 

Proof: We need to be careful near the boundary at zero. Begin by extending onto R, for example 
by setting (j)(x):= \$" (0) x 2 + 0' (0) sc for x < 0. Then G C 2 (R) so we can apply Ito's formula to 
X. ) which gives that for almost all u) G 



M + 



1 



+<TI 



f'xDdwl 



for all i G [0, T]. We can replace t by t Ar, noting that ^X t 1 Ar J = (f> (X tAr ), since (f)(0) = (f) yX^: 
= 0(9), and 



t pt/\T rt 

s 

./0 J!Ar 



' ' V, 1 ) <!s - I o' (Xl)ds+ j <p' (8) ds = j ' I <• 1 



X l s ) ds 



and likewise for the other integral terms, therefore 

# (X, 1 ) = v4l (X\) + l - (a 2 M + aj) 0" {X}) dt + a M <P' {X}) dW t + orf {X l t ) dW x t 



and the result follows by taking expectation over Q 1 and noting that E 1 [(j) {X})] 



\vt) 



□ 



Remark 2.4 (Connection with filtering) If we treat X 1 as a signal process and W an obser- 
vation process, then the Zakai equation (see for example pQ) for the conditional distribution of the 
signal given the observations is exactly ([!]). 
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2.2 Properties 

For a fixed uj G Vt the process X 1 which is not killed upon hitting the origin is simply a iV-Brownian 
motion started at Xq with drift t h-> \xt + auWt- We can therefore explicitly write the f2 -density 
of X 1 as: 

ol / \ tjI i \ t \ 1 j (x-x - fit- a M W t f \ 
P t (x; x ) = P t (x; x ) (u) = exp I \ . 

It is trivial to see that if X\ G S for a time t G [0, T] and a measurable subset S C (0, oo), then one 
must have X\ G 5. This simple observation allows us to show that v has a density process: 

Proposition 2.5 (Existence of the density) There exists a full subset of O on which for every 
t G [0, T], Vt has a density process Vt G L p (0,oo) for all 1 <p < oo 7 urei/i ||Vf|| < ||Vq|| . 

Proof: From the above observation, there exists a full subset f2 on which 

v t (S) < F 1 (Xl eS)= J P/ (x; x ) V (dx ) dx 

for all i G [0, T] and S C (0, oo) measurable. Therefore if S Q (0, oo) has zero Lebesgue measure 
then vt (S) = 0. It follows by the Radon-Nikodym Theorem that vt has a density for every t G [0, T], 
and that this density Vt satisfies 

rod 

V t (x) < / P* {x;x )Vo(x )dx (3) 
J o 

for almost all x G (0, oo). 

Bounding Vq by HV&lloo in we get V t G L°° (0, oo) with WVW^ < WVqW^. Suppose 1 < p < oo, 
then Holder's inequality applied to ([3j gives 

/■oo 

n(^) P < / P/(x;x )T/o(^o) P ^ . 

J 

Integrating over x G (0, oo) and noting that J °° P t (x; xo) dx < 1 gives ||Vt|| p < || Vb|| p . 

The most important feature of the probabilistic solution is the following estimate: 

Theorem 2.6 (The (3 + /3)-Condition) Let V$ G L 2 (0,oo) be bounded and a be defined as in |I|). 
Then for every (3 G (0,n/a — 1) C (0, 1) there exists constants B > and 7 G (0, 1) such that for 
all t G [0, T] and e > 



□ 



E 



< -e 3+ ^. 
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A proof is presented in section 3. 

In Lemma 3.5 of jl] this result is presented for an initial density supported on a finite interval 
away from the boundary. This assumption allows the authors to drop the singular i~ 7 factor. 
However in some respects this is not a natural initial condition. For any t > 0, Vt takes positive 
values on all of (0,oo), so if we stop the process at time to and then restart from the initial density 



Vf we still have the estimate on E 



for t > to j despite Vt not being supported on a finite 
setting we can expect a singular time factor to appear 
> (Theorem 3.2 of [9j). Therefore positive values of 



interval away from zero. In this more genera 
since for every time t > 0, Vt (zc) — > as x ■ 
Vq close to zero must decay instantaneously as the system evolves in time. Since 7 G (0, 1), i~ 7 is 
integrable over [0, T] and so this factor does not present technical difficulties in later proofs involving 
L 2 -integrability. 

{XlXl). Under 



2.2 



we prove Theorem 2.6 by considering the process 



Using Proposition 

PxP 1 this process is a two-dimensional Brownian motion that is killed upon exiting the quadrant 
(0, oo) 2 and has components with correlation 

CM 



a 



M 



With a as in l\2h, the map F a 



I 2 defined by 



Fc, : x )->■ 







1 



X 



(4) 



transforms X to a Brownian motion with uncorrelated components that is killed upon exiting the 
wedge 

(0,oo) 2 = {(r,6») : < r < 00, < 9 < a}. 

We then make use of explicit formulae in [5], |12| to estimate the hitting probability of a small 
neighborhood of the apex of the wedge. The kernel smoothing technique introduced in |3] and 
adapted from [11] then allows us to relate this hitting probability to the behaviour of the solution 
near the absorbing boundary. 



3 Proof of Theorem 12.6 

Let us write ¥ xliX2 to indicate that Xq = x\ and Xq = x%, that is we start the two-dimensional 
Brownian motion X = (X^,X 2 ) from {x\,X2). From Proposition 



2.2 



we know 



[vt{S x )v t {S 2 ))=^ 



Xl,X 2 



E 



L xles l7 xfes 2 



(5) 
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for all S±, S*2 C (0, oo) measurable. A change of measure allows us to eliminate the drift term in the 
dynamics of X: 

Lemma 3.1 (Eliminating drift) Let P be defined on (£l,Tt) by the Girsanov transformation 



dF 7 ( \ 



so that under P x P 1 X is a two-dimensional (fi x Q^-Brownian motion with zero drift. Let E 
denote the expectation under P. // there exists a constant C > such that 



E 



MM 



< 



fir /2a 



2+w/a 



for all t G (0, T] and e > 0, i/ien Theorem 2.6 holds. 

Proof: Following the method in Lemma 3.5 of [4J, for any Holder conjugates 1 < a, b < oo we have 



E 



^(0,e) 



E 



Z t (ij)-v t (p,ey 



< t[Z t (ii) a ]»t [v t (0,e) 

< E[Z t (^) a ]"E U(0,e) 



\26 



< E[Z t (^) a ]° 



£7r/2a! 



2+vr/a 



E 



a-l 



1 ' C -.,!«- 



W 2a 



2+tt/q 



where one should note that v% (0, e) < 1. Hence 

(a -1) (a -2)fi 2 r 



E 



< 



exp 



2a 2 M a 



C 



£ 2+ir/a[ ; 



therefore we can choose 6 sufficiently close to 1 so that the result holds since a E (n/2,ir). 



□ 



Remark 3.2 Notice that we cannot take 6=1 above. This change of measure technique eliminates 
the possibility of taking /3 = ir/a — 1, except in the case when fi = 0. Although this may seem to 



weaken the potential regularity available near the default boundary, in proving Theorem 1 . 1 we make 



use of Lemma 4.8 which requires j3 to be strictly smaller than this maximal value. It is therefore 



unclear as to the precise behaviour of V at the critical value /3 = n/a — 1. 
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The remainder of Section 3 is devoted to showing the hypothesis of Lemma 3.1 holds. Clearly it is 
sufficient to take fi = and P = P. 

Conditioning on the start-point of X: 



E 



^(0,e) 



o jo 



E 



Xl,X2 



E 1 



L x t 1 e(o,e),x t 2 e(o,£) 



Vq (xi) V (x 2 ) dxidx 



2- 



Now the region (0, e) is mapped into a region of the form {(r, 9) : < r < ce, < 6 < a}, for some 
numerical constant c, under the transformation F a . As the estimate we are demonstrating is un- 
changed under the map e i— > ce (modulo multiplicative constants), it is no loss of generality to take 
c = 1. Using the fact that Vq is bounded, that the Jacobian of the transform F a : M? — > M? is a 
numerical constant depending only on p, and formula (8) of [5] (corroborated in |12j). we arrive at 

tittOq 



E 



f t (0,e) 



oo pa ret re 



CI 111 —e- — 
lo Jo Jo Jo t 



E 

n=l 



mrb 
sin sin 

a 



a 



'-Inn /a i^f) drdBdr^ 



o 



C 



e ir/a+1 



4 E 



+ 1 



f \ m+mt/ct 

— ) dr 
2tJ 



for C > a numerical constant. As sup x>0 {x c e x } = c c e c we have 



E 



e w/a+1 

1 n,m=0 



(m + + l) (m + ^f) m+ ^ e-( m+ ^) 

(2n + l) 2 m\F (m + + 1 



(6) 



So we are finished provided the double sum on the right-hand is finite. Using Stirling's approxima- 
tion, the inequality 

r(q + c+l) r(a+l) 

r(6 + c + i) - r(6 + i) 

for a,b,c > and the asymptotic behaviour T (m + -k /2a + 1) ~ m\m n ' 2a for large to, gives that 
the summand in ^ is dominated by a constant multiple of 



(2n + If m w /a«+i/2 



for large m and m + n. Since 7r/2a + 1/2 > 1, this expression is summable over m and n and so we 
have the result. 



□ 



10 



4 Proof of Theorem 



1.1 



We begin by proving the regularity result for the probabilistic solution. Let us write ps for the 
transition density of a Brownian motion in R with zero drift and unit volatility, that is 



for z G R, 5 > 0. 

Definition 4.1 (Heat kernels) Let 5 > 0. T/ie absorbing and reflected heat kernels are defined 
to be 

Gs 0, y) := (a? - y) - pa (x + y) 
G l{x,y) '■= P&{x-y)+Ps{x + y) 

for x, y > 7 respectively. The smoothed density and reflected smoothed density are then be defined 



for t £ [0, T] and x > 0, respectively. 

Now take any smooth compactly supported function <f> G C|? (0, oo). Then for every 5 > 0, 
x I—)- Ts4>(x) (that is the function obtained by convolving the heat kernel with (j)) is in C test , so 
putting Tg(j) and into ([T]) we arrive at the following strong equation for the evolution of the 
smoothed density: 

Proposition 4.2 (Strong smoothed evolution equation) For 5 > 0, t G [0, T] and x > 




(7) 



to be 





T ,s (x) + 



'1t^(x)1 ds-ciM f t T' s>s (x)dW s 







(8) 




where R t $ is the remainder function: 
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Proof: This is simply a matter of observing that 

(<t>',T t , s ) = -(*, [n )S y) 

One also needs to apply the stochastic Fubini theorem [3] to switch integration in the space dimension 
with integration with respect to W (and likewise with respect to time). To this end it is enough to 
note that the tails of Tfg and TJ S and their derivatives decay exponentially, hence the conditions of 
the (stochastic) Fubini theorem are indeed met. 

□ 

Differentiating equation (|8j n times (or integrating once for the case n = — 1, note the result is 
no different) and applying Ito's formula to \ T^j gives: 

+4/ 2 * + ^R^dt (9) 

-2a M T$T% +1) dW t + 2a M T^R^ 1) dW t . 

From this equation we shall proceed to the proof by induction in section 4.2. Before we continue we 
introduce several technical lemmas. 



4.1 Technical Lemmas 

Lemma 4.3 There exists a full subset of Q on which x~ c Tt^$ (x) — > as x — > for every c £ (0, 1), 
S > and t € (0,T]. 

Proof: This is straightforward since for x < 1 we have 



< T t>5 (x) 



Jo 



exp 



25 



exp 



(x + yY 
25 



Vt (y) dy 



6V2tt6 Jo 
D s x 



y exp 2S~\dv 



for a numerical constant D$ > 0, and where we have used the estimate 1 — e z < z. 



□ 



12 



Lemma 4.4 Let n > be an integer. Suppose that for k = 0, 1, . . . , n 



lim inf E 

5->o 



o 



dt 



< oo. 



Then for almost all (u,t) G x [0,T], the density Vt is n-times weakly differentiable with 

rT 



E 



o 



(fc) 



< oo 



(10) 



for k = 0, 1, . . . ,n. 



Proof: Let (•, •) denote the usual L 2 (0, oo) inner product: (/, 5) = / °° / (x) g (x) dx . Let {(j)i} c *L l 
be a basis of L 2 (0, 00) such that each <f>i is smooth and bounded. 

Firstly for any <fi G C°° (0, 00) and m > k > 1, lim^o (%-^/2^ (x)) = 0, since d^w m _p/2 is a 
linear combination of w i _p/ 2 ~~ ^ 0, for m — k < i < m. So for > 1, integrating by parts gives 



Applying Leibniz's rule to (w k _i_pi 2 <j)) , we see that 



(-!) fc \(/>,™k-l-p/2 T $) ~ 1 j™{ T *'' 5 ( x ) (^fc-i-/3/2<A) (fc 1} } 



(fc-l) 



Ijm \r t>s (x) (w fc _i_/3/2^) (fc 1} (^)} = Cjfclkn{T ti i ( x ) w -P/2 (x) <f>(x)} 



for some numerical constant C k , and the right-hand-side vanishes by Lemma 4.3 since w_p/ 2 (x) 
O (x - ^/ 2 ) for small x. Hence we conclude 

( {w k _ x _ p / 2 (j)) (fc) ,T t>s j = (-l) k (fawk^p/zltj) ■ 

The remaining case k = is trivial. 

From these results and Fatou's Lemma we have 



E 



OO 

E 

i=l 



>fc_l_^/20i) (fc) ,^ 



eft 



E 



_y 00 

J0 i=l 

oo 

E 



< lim inf E 

5-5-0 



= lim inf E 

5->o 



dt 



T oo 



eft 



I T^( fc ) 



< 00 
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for k = 0, 1, . . . n, therefore we can define the processes 

oo 

W tJk (x) := J2 ( K-1-/J/2&) W , V t ) 4>i ix) 

which satisfy 



8=1 



E 



\\Wt,k\\Ut 



< oo. 



For k = 0, let /q = W-i-p/2 and take any ^ G C|? (0, oo), then we calculate 



ip, -rWt,o 
Jo 



i=l 



and so Vt = Wto/ fo- Hence we have demonstrated (10) for k = 0. 



Take k > 0, and set /& = 1-/3/2 ■ Using Vt = W^o/Zo we can compute 



Jk 



i=l 



1 



E ((/^) (fe) ,|)(/o^,14) 



V' 



, JO 



fo 



*.J=1 
00 

E ((/*&) 
i,i,m=l 
00 

j j',m=l 

E (-i) fc (v,4 fc) )(^,0m)(/o^,F t ) 

j,m=l 



Jk 



hence (— l) fc Wt^/fk is the fc th weak derivative of 14, which completes the proof. 



Lemma 4.5 There exists a full subset of O on which 



T (n) 



n 



w 



(n) 



n-l/2-p/2 T t ',8 



and 



14 



w 



n-fl/2- L t, 



T$>(x)T% +2 \x)dx = 2 



w 



rp(n) 

n-/3/2 I t ,8 



It' 



-4 n 



it: 



n-/3/2- L t,6 
» 



n-l/2-/3/2^ t ,J 



+2(n-- n- — 



tv 



71-1-/3/2-*+ J 



/or every n > —1, (5 > anc? i G [0, T\. 

Proof: The result is an integration-by-parts exercise provided the following limits vanish: 

So { Wn ' 1 /2-/3/2 0*0 r +. ( ,? 0*0 } 



Iim{^ /2 (.)r5)(x)r5 +1 )(x)} 



for n > — 1. If n > 1 then u> n _i/2-,3/2i w n-/3/2 — ?• as x — )• 0, and the result is immediate. In the 
case n = 0, Lemma |4.3| gives (x) = O (x) which is sufficient for the limit to be zero. This also 
implies T^g ^ (x) = O (x 2 ), so the result holds for n = — 1 too. 



□ 



Lemma 4.6 With ps defined as in |?|) ; for every integer n > there exist constants cf for n/2 < 
i < n such that for every 5 > and x > 



d n ps 
dx n 



£ c^-^psix). 



n/2<i<n 

Proof: This is a simple inductive argument. 

Lemma 4.7 For n = —1, 0,1, ... ,N and (3 < ir/a — 1 



□ 



lim inf 

<5->0 



(n) 



< oo. 



Proof: The cases n = and n = — 1 are simple. Firstly observe from Definition 
HVblL- Therefore for all 5 > 

|h-0/2Zb,*||2 < IIVbHoo / X-Pe~ 2x dx < OO, 



4.1 



that llTr 



o,<5 1 



< 



15 



and so we have the result for n = 0. The case n = — 1 follows by noting that < T^ 1 (x) < ((VoH^ x 
and applying the same argument. 

Now fix n > 1. We begin by noting that d™G$ (x,y) = dyG r ^ (x,y), where G r " = G if n is 
even and G""™ = G r if n is odd. Splitting the range of integration and integrating by parts in the 
definition of T $ gives 



2 1 /"OO 

£ ^^-Gf (x, x/2) V^ l) (x/2) + / Gf (z, y) F (n) (y) (11) 



+ 



x/2 



d n x G s (x,y)V (y)dy. 



We consider the weighted L 2 -norm of these three components. 



Firstly by Lemma 4.6 there exists a numerical constant D > such that 



■yi 1 



Gf (x,x/2)^ « (x/2)w n . p/2 (a 



< D 



,2i-n+l+i r i-l/2 e -| 3 



V W (x/2) 



■i-j8/2 



X X 



-<i<n-l-i 



therefore integrating over x £ (0, oo) and applying the Cauchy-Schwarz inequality 

2 



d^Gf -/2) F (i) (-/2) u^/j, 



X 







<2D ^ 

< D' ^ ^/2 2 2i-/3 



4i+2 r 2,-l e -| Jdx 



y (i) (-/2)^ /2 



^0 Wi-p/2 



(12) 



as (5 — > 0, where D' > is a further numerical constant. 

For the second term we apply Cauchy-Schwarz to obtain 



x/2 



GT(;y)vi n) (y)dy 



< 



Gf (-, y) w n _ m (y) 2 Vt ] (y) 2 dy\ll Gf (•, y) y-^e^dy 



x/2 



< 



4 . 2 2n-P x -2n +P U ps {x _ y) {y f y(") {y f dy \ J f p g ( x - y) ^ d y . 



x/2 
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For 5 < 1 we have f^PS ( x — y) e2y dy < D'e 2x , for some numerical constant D' > 0, therefore 



I G r 5 n (;y)vi n) (y)dy 

/2 



< 2 2n-p + 2 D , > I p s ( x -y)dxw n ^ /2 (yfvi n) (y) 2 dy 
Jo Jo 

2 



< 2 



2n- / 3+2£ ) ' 



(n) 



(13) 



Finally using the fact that V$ is bounded and Lemma 4.6 there exists a numerical constant 
D" > such that 



x/2 



d^G s (x,y)V (y)dy 



n/2<i<n 



and therefore 



•/2 



™n-0/2 / (-,!/) Vb (y) dy 





< 2{D") 2 \\V \t E / ^" 2n+2 ^ 1 ^^dx 

2 n/2<i<n 

< £,"',$(1-40/2 (14) 



as <5 — )• 0, since /3 < 1. Here JD"' > is a further numerical constant. 



Results (12), (13) and (14) give the required result by (11). 



□ 



Lemma 4.8 For n > and j3 < ir/a — 1 



lim inf E 

<5-s>0 



(n+1) 



W n -P/2R S ,& 



ds 



0. 



Proof: Interchanging differentiation and integration and applying Lemma 4.6 gives 



/■OO 2 2 

< D I (x + y) 2l - n - 1 S- l - 1 / 2 e- X 26e- V ^V s (y)dy, 



=±i<i<n+l 
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with D > a numerical constant. Let r/ £ (0, 1) and split the region of integration at y = d^ 1 ' n \ 

y 1 y 2 l 

then for 5 < 1 we have 5~ l < 5~ n ~ l and e~zs < e~~ eT l l^ on y > 5^^^ so that 



Wn-p/2 (a?) 1) (x) 



< 



=±i<j<n+l 



5<i— n)/a 



/ . \2i-n-l -2? -s£ 
(x + y) e 2 e 4 



iVblL^-^e-V^^ (x) (is) 



^<i<n+l 



2%-n-\ 



The first term on the right-hand side of (15) vanishes as 5 — > 0, hence we have 



lim inf E 

5-5-0 



(n+l) 



ds 



< lim inf D' ^ if (x + 1) 

x s (l- V )(2i-n-l) 6 -2i-l 5 ^(l- V )(3+l3+e) 6 l/2 § n-l3/2^ (jg) 



where we have used Theorem 2.6 with /3 replaced by /3 + e < ir/a — 1 for small e > 0, and D' > is 
a numerical constant. (Note as we integrate over s E (0, t), the singular factor of s~ 7 from Theorem 



2.6 integrates to a finite value because 7 G (0, 1).) The 5 term in (16) is bounded by 

5 \{l- V )e-n(2n+b/2+l3/2) 

for small 5. We choose 77 small enough so that 

J(l-r ? )e>? ? (2n + 5/2 + /3/2) 



then the limit on the right-hand side of (16) vanishes and gives the required result. 



□ 



4.2 The inductive proof for the regularity result 
Initial case: 

The initial case is when n = — 1, and the result we shall demonstrate is: 

-T 



lim inf E 

(5-s-O 







w_2-p/2T; s 



(-1) 



dt 



< 00, for every (3 < n/a — 1. 



(17) 



By truncating the range of integration we have 
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Tt s 1] (x) = / / G s {z,y)dzV t {y)dy = v t (0,x) + I t>5 (x) 
Jo Jo 

where I ti s (x) := J Gs (z, y) dzVt (y) dy. From Theorem 



2.6 



it is clear that 



E 



Uo 



I 1 1 z 

\W-2-l3/2^t(0,-)\\ 2 dt 



< oo 



and therefore we are done provided lim inf K Jq \\ w -2-i3/2dt,s\\'^ dt < oo. Using the inequality 



G5(Z ' y) -^& 6XP 
we have a numerical constant C > such that 



{z-yf/25}, 



hs (x) < C5- 3/2 x 2 e- x2 / 25 / (y + 2x) e - y2/25 V t (y + 2x) dy. 

Jo 



and so 



i2 
I2 

2 roo poo roc 



c 



< 



5 3 



roo roo roo x^+yj+yj 

x V (yi + 2x) (2/2 + 2x) e 2* 7 t (2/1 + 2a;) Vt (y 2 + 2x) d yi dy 2 dx. 
Jo Jo Jo 



Now let 7] E (0, 1). Splitting the region of integration on (0, c^ 1-1 ?)/ 2 ) 3 C (0, oo) 3 and its complement, 
we see that on the complement the triple integral on the right-hand side above is bounded by 
a multiple of e -1 /^ , and so vanishes as 5 — > 0. Therefore we need only consider the region of 
integration (0, S^~ v ^ 2 ) 3 , and so we have a numerical constant G > such that 



lim inf E 

5-5-0 



< C'liminf^ 3 ^)/ 2 / E 
<S-H) 



dt 



where the second line follows by Theorem 2.6 and e > is sufficiently small. Clearly we may choose 
77 > small enough so that the exponent of 5 is positive, whereby the right-hand side vanishes and 



we have (17). 
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Inductive step: 

Assume that for some — 1 < n < N — lwe have 

r-T 



lim inf E 







rr(k) 
W k -l-f3/ 2 T t \ s ' 



dt 



< oo, for every /3 < n/a — 1 



(18) 



for k = — 1, 0, 1, . . . n. (We demonstrated the case n = — 1 above.) The proof will be complete if 
we can show condition ( |18[ ) holds for k = n + 1, since then by induction the condition holds for 
all fc = 0, . . . , iV and we can apply Lemma 4.4 to give the regularity result of Theorem |1.1| for the 
probabilistic solution v. 

Let us begin by multiplying throughout equation ([9| by the square of the weighting function 

and arrive at 



•iv 



n-p/2 



and integrating over (0, oo). To do so we apply Lemma 



4.5 



o-fE 



w n -0/2T SiS 



O+l) 



rn{n) 



E 





n 

V 2 



E 



-2/dE 

+2 (a| f + <r?) E 
-4 (<4 + a? 



t 



(n) 



ds 

2 



(19) 



o 



^-1/2-/3/2^5 



(n) 



c/s 







(n) 



n ~ 2 ' E 



t 



ds 







(n) 



+2 a^ + a? n 



/3 



n 



+<r^E 
+2/iE 



t roo 



(n+1) 



/3+1 



E 



(n) 



c/s 



JO 



...2 r Wn(i+l). 
w n-/3/2 1 s,5 U s,5 as 



Eliminating the negative terms from the right-hand side and applying the Cauchy-Schwarz inequality 
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gives 



a 2 E 



o 



< 



2 



(n+1) 



™n-/3/2-L ,6 
t 



-2/xE 



o 



rr{n) 

w n _p /2 T y s J 



ds 



ds 



-2// (n-fl E 



+2 (a| f + af ) E 
+2 (4, + of) 



(n) 



2 



ds 



w 



n-(3/2J- s ,& 



ds 



n 







n 



2 

P+l 



E 



rr(n) 



(is 



+a 2 M E 



+2|a*I i E 



n-f3/2R s ,6 



(n+1) 



(is 





m(n) 


2 










ds 


}"'{« 


[/' 




.J 


2 




yo 





(n+1) 



1/2 



(20) 



We are done provided the seven terms on the right-hand side of inequality (20) have a finite limit- 
infimum as 5 J. 0. This is true of the second to fifth terms by the inductive hypothesis (18), and with 
this the final two terms vanish by Lemma |4,8| From Lemma |4,7| the first term remains finite in the 
limit-infimum, hence we have shown the probabilistic solution v satisfies the regularity properties in 
Theorem 11,11 



4.3 Uniqueness 

It remains to show that the probabilistic solution is indeed the unique finite-measure valued solution 
of ([T]). Suppose we have another finite-measure valued solution v. Define the signed-measure valued 
process u as the difference between the solutions: 



u t (S) := v t (S) - v t (S) . 

By linearity of ([T]), provided that the technical lemmas in Section 5.1 hold if we replace v by u, then 
we arrive at (|19|) and the same limit results hold. Whence in the case n = and /3 = we have 



<r?E 



\ U 's,s\\2 ds 



-E 



+ 2 (a 2 M + &j-n) I E 
'o 



\Us,6 



ds + f(t,5), (21) 



where Ut s denotes the function obtained by applying the kernel smoothing operator to m and 
/ (t, 5) —> as 5 —7- uniformly in t £ [0,T]. Observe that the initial condition vanishes because 
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uq = 0. Since the left-hand side of (21) is non-negative, we have 



lim inf E 

5->o 



o 



4.4 



is 



for all t £ [0, T] by Gronwall's Lemma. This is then enough to show u = provided Lemma 
applicable for u. 

It therefore remains to show the working in Section 4.1 is valid for v replaced by u and the 
specific choice N = 0, /3 = 0. Firstly observe that by taking expectation over ([!]), the (determin- 
istic) measure- valued process (E [i*t]) t6 [o t] satisfies the (weak) heat equation in the half-line with 
zero boundary condition and zero initial condition. As ut is a finite-measure, the uniqueness of 
solutions to the heat equation forces E [ut] = 0, and so E [i/J = E[z^]. Likewise, by considering 
the two-dimensional heat equation (with a correlation term), we conclude that E \yt (Si) vt (£2)] = 
E [i?t (Si) i>t (£2)], f° r a H measurable Si, £2 C (0, 00). Therefore by the Hahn Decomposition Theo- 
rem we can write 



E 



u t (0,e) 



< E 



+ E 



2E 



This is enough to guarantee Lemma 4.8 holds for u. As Lemma 4.3 may no longer be true for u, the 
remaining facts that need to be checked are whether the limit lim x ^o Ut,s (x) = from Lemma 
and the limits 



4.4 



YmxU t fi (x) io_i/ 2 (x) = 
limU t;S (x)Ul >s (x) = 



from Lemma 4.5 remain valid. However it is easy to see that 



\U t ,s(x)\ < / {p s (x - y) - ps (x + y)} \ut\ (dy) 
Jo 

\2 "l / 



< 



1 



y/2n6 Jo 
2x 

5V2nS jo 



(x-vY 

exp <^ - 
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1 - exp <{ — S 



\ut \ (dy) 



00 J (x-yf \ 

ex p \ — } y M (dy) 
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where \ut\ < vt + ^t- For x < 1 the integrand is bounded and hence the integral is finite since \ut\ is 
a finite measure. We therefore conclude that Ut t s ( x ) = O (x) as x — > 0, which is enough to complete 
the proof. 



□ 
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5 Proof of Theorem 



1.2 



This theorem relies on the fact that the estimate in Theorem 2.6 is sharp. Specifically, Lemma 5.1 
of [6] states that if Vo is supported on a finite interval away from zero, then for every r] > there 
exist C (rj), Eq (rj) > such that for all < e < £q (v) an d t £ [0, T] 



E 



M0,e) 2 l > C(r])e 2+7T/a+71 . 



(22) 



It is clear that if Vo < Vo then the empirical measure v arising from starting the system from Vo 
satisfies v% < Vf Since we can certainly find a constant a > such that Vo := Vol(a )0 o) i s non-zero, 
we see that (22) remains valid if we drop the assumption on the support of Vq. 



For a contradiction, suppose that we have /? > ir/a — 1 and a k £ {0, 1, . . . , N} for which 

rT 



E 



Then in particular 



E 



o 



T rl 



(k) 



dt 



< oo. 



^0 



■ x 2(k-l)-{3y(k) {x) 2 dxdt 



< OO 



for any (3 < f3, so let us fix /J £ (tt/cx — 1,1). Also take e > small enough so that we have 
/3 - e G (vr/a - 1, 1). Then for fc > 2 



c 2( fe -2)-^ +£da . 



and since 



(v-^^-^a)) 2 = (jy t {k) (y)dy 



< f 1 x W-l)-Py(k) {x) 2 dx 

Jo 

< ! x- 

Jo 



Xk-1)-Py(k) {x) 2 dx 



y 

^2fc+3 



2k -3- j3 



(23) 



applying the triangle inequality, multiplying by x 2k 4 13+6 and integrating over Q, x [0, T] x (0, 1) 
gives 



E 



T rl 



^0 



(x) dxdt 



< oo. 



Repeating this argument eventually yields 



E 



f [ x~ p+e Vl (xf dxdt 
Jo Jo 



< oo 
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for all allowable e. If we return to (|23j) with k = 1 and take the upper limit of the integral to be x 
and the lower limit 0, then since Vt (0) = we arrive at 



E 



T rl 



JO 



- 2 -P +£ V t (x? dxdt 



< oo. 



Repeating this manipulation a final time gives the following for x < 1: 



E 



v% (0, x) dt 



< E 



T /■! 



J 



x -2-/3+ E ^ ^2 dxdt 



y 2+p - e dy 



J J 



D(P,e) x 



with D £ (0, oo) a constant. This is a contradiction of (22) for the case r/ = (1 —n/a) + 

/3 — > 0. Hence we have Theorem j 1 . 2 1 

□ 
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